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L INTRODUCTION 

The goal of this research is to gain quantitative fundamental 
understanding of the properties and behavior of clusters of metal atoms 
as a function of the number of atoms in the cluster. Among the properties 
of interest are the various possible equilibrium structures a cluster of 
atoms may adopt as each atom is added, the relative stability and degree 
of order of these equilibrium structures, their characteristic fundamental 
modes of vibration, and the magnitude of the energy barriers regulating 
the rate of interconversion of one equilibrium cluster structure to another 
one. 


II. SCIENTIFIC PROGRESS 


During this reporting period, scientific progress was made in 
two different areas. Continuing progress was made on developing the capa- 
bility to calculate gradients of the matrix elements of the effective core 
potential. Particular emphasis has been given towards finding a suitable 
algorithm for the efficient processing of the derivative integrals after 
use has been made of the combined translational and rotational invariance 
properties of the integrals to reduce the integral calculation to just the 
truly linearly independent derivative integrals. 

The most significant progress during this period, however, con- 
sists of the development of an entirely new approach to the calculation of 
electronic potential energy surfaces. The details of this new approach are 


described in Appendix A in the form of a preprint of a manuscript that 
has been submitted for publication to The Journal of Chemical Physics, 

This new approach is based upon the atom-superposition-and-’electron" 
delocalization (ASED) model of chemical binding. It has the advantage of 
dividing the electronic energy int*;^ physically transparent groups of terms. 
Moreover, whereas conventional methods approach a molecular calculation 
without recourse to the information already calculated from the component 
atoms, the structure of the ASED energy lends itself to the fullest utili- 
zation of the properties already available about the atomic solution in order 
to simplify the calculation of the molecular energy. This new approach has 
prospects of reducing even by an order of magnitude the effort required in 
the calculation of the molecular energy, and thereby constitutes a major 
step towards accomplishing the goals of this research explained in the 
introduction. 
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A new approach to the calculation of potential 
energy surfaces 
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ABSTRACT 

A new decomposition of the molecular energy is presented that is 
motivated by the atom-superposition-and-electron-delocalization physical 
model of chemical binding. The energy appears in physically transparent 
form consisting of a classical electrostatic interaction, a zero-order 
two-electron exchange interaction, a relaxation energy, and the atomic 
energies. Detailed formulae are derived in zero- and first-order of ap- 
proximation. The formulation extends beyond first order to any chosen 
level of approximation leading, in principle, to the exact energy. The 
structure of this energy decomposition lends itself to the fullest util- 
ization of the solutions to the atomic sub-problems to simplify the cal- 
culation of the molecular energy. If non-linear relaxation effects remain 
minor, the molecular energy calculation requires at most the calculation 
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of two-center two-electron Integrals. This scheme thus affords the pros- 
pects of substantially reducing the computational effort required for the 
calculation of molecular energies. 
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I. INTRODUCTION 

In moleculan problems involving large numbers of electrons* as in 
the case of molecules consisting of many transition-metal atoms and many 
carbonyl groups, the application of the standard first-principles methods 
of molecular quantum mechanics becomes Impracticable, In these cases 
one must resort to more approximate methods^^*^^ that, while practical, 
retain the prospects of yielding a realistic physical description of 
chemical binding. 

One semi-empirical method that has been actively applied to large 
electronic-structure problems is the atom-superposition-and-electron- 
delocalization (ASED) method yf Anderson^®^ This method has the desir- 
able property that, in addition to the molecular one-electron energy levels, 
it also yields an approximate molecular total energy as a function of the 
function of the positions of the constituent atoms. The application of 
this method to large molecules has been shown to yield useful predictions 
of molecular structure, force-constants, and relative bond strengths 
The application of the ASED method ranges from the study of the structures 
of clusters of transiti on-metal atoms to the study of molecules chemi- 
sorbed on metal-cluster surfaces^®^, to the study of organometaVlic com- 
plexes(9'10). 

The ASED method is a semi-empirical method. The molecular orbital 
energies are obtained by what is in essence the extended- Htickel method^®*^^^ 
The ASED molecular energy MUatz is obtained by combining an energy 
that is derived from an apnroximation to the integral Hellmann-Feynman 
force formula^^’^^^ with the sum of molecular orbital energies (the extended 


Hlickel energy). 
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The validity of the ASED wethod rests mainly on the evidence provided 
by the numerous useful applications that have resulted from its use. There 
exists, however, no first-principles derivation of the ASED method that 
would show how this method fits as a particular step into a well-defined 
hierarchy of approximations to the exact molecular energy. The absence of 
such a theoretical basis has been the source of skepticism about the validity 
of the ASED method. Recently, Anderson has presented work aimed at elucidating 
the theoretical basis of the ASED method. While the physical reasoning 
in this work, based on a time-dependent approach to chemical bonding, is 
compelling, the formulation of the analytical aspects of the method, espe- 
cially for the many-electron cases, remain heuristic. 

In this paper we present a first-principles decomposition scheme of 
the molecular total energy using the conventional procedures of molecular 
quantum mechanics. We have attempted to adhere, as closely as we found 
possible, to the ASED physical model of the stages in the bonding process 
to set up the basic steps from which then flow, as a consequence, the 
definitions of the terms in the decomposition of the energy. This energy 
decomposition is formulated so that it may, in principle, be taken to any 
chosen level of approximation leading up to the exact energy. There are 
some inherent differences in a semi-empirical approach and a first- 
principles approach such as the present one. Without attempting to resolve 
this issue, we find that in zero-order of approximation the present energy 
decomposition does indeed have an analogy in the form of its terms to the 
semi -empirical ASED energy aMiUz, This analogy in form, however, does 
not extend beyond the zero-order level. The inclusion of the first-order 
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correction terms yields a different energy formula that, nevertheless, 
corresponds physIcaTly even more closely to one of the main effects 
envisioned 1n the ASED physical picture of the chemical bonding process, 
namely the relaxation effects subsequent to the ''Igid-atom superposition 
stage. The main themes of the ASED physical model of chemical bonding 
are found to reemerge In generalized form even on Including higher order 
corrections to the energy decomposition. 

II. THE ASED SEMI-EMPIRICAL METHOD 


Procedurally, the ASED method consists of two steps. The first 
step Is the calculation of pair-wise Interaction energies among the atoms 
in the molecule using the formula 



The summations range over all unique pairs of atom labels A and B, and 
the symbols Z^, and for example, refer to the position 

vector, nuclear charge, and atomic charge density of the A'th atom re- 
spectively. The spatial coordinate vector Is written as The second 
step is the calculation of the molecular orbital energies, {e|^, k *= 1, 
...» n), by solving the secular equation 

^ L (^^v,A;v,B “ 

A PCS} 

A 

vCSg} k “ 1 , . . . , n N (2) 

The restrictions and approximations to this equation are the same as in 
the extended Huckel method. ^ The basis set on each atom, 5 is 


u,A 




V.A;k 


S 0 
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restricted to a uiinitnal set of valence atomic orbitals. The 

are the elements of the overlap matrix, The hamiltonian matrix elements 

are defined as 


^p,A;v,A ' " ^y,A 


(3a) 


The I A is the ionization potential for removal of an electron from the 

]Xfn 

y'th orbital of the A'th atom. The K factor, based on experience with 
first-row diatomic molecules, is defined by Anderson' * as 

K - 2.25 . (4) 

where D is the distance between the centers. Note that we shall be 

using Hartree atomic units in all the equations. When required, the 

0 

following conversion factors are used; 1 bohr * 0.529177A, 1 hartree = 
27.21165 e.V. The approximate total energy of the ASED method is ob- 

.. /ii) 

tained by combining the extended Huckel energy' ' 

^EH " S k ^k ’ 

k=1 

where (n_„)^ is the occupation-number of tlie k'th molecular orbital, 

OCC K 

with the Ep interaction energy, Eq. (1) as follows 

E=ER.EgH (6) 

The physical picture that guides the formulation of the ASED 
method envisions chemical binding as the synthesis of two distinct 
processes. The molecule is envisioned as being built-up, first, 
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by the process of bringing together "rigid" and "non-exchanging" atoms, 

The electrons of each atom are not allowed, at first, to be perturbed by 
the proximity of other atoms, and neither are the electrons allowed to 
exchange between the atoms. Once the atoms are thus positioned, the 
second stage of the molecule-building process is envisioned as allowing 
the electrons to "delocalize" from the parent atom to all other atoms. 

In order to isolate the contributions to the energy of the exchange 
process, effects on the electrons due to fields from neighboring atoms, 
according to Anderson, are to be neglected during this second stage, 

The first process in the build-up of the molecule Is envisioned as being 
energetically unfavorable for bindingi the second process is envisioned 
as being energetically favorable. The balance of these two is envisioned 
as leading to the equilibrium configuration of the atoms in a molecule. 

It is convenient to separate the issue of the validity of the energy 
fliuatc used in the ASED method from that of the usefulness of the physical 
picture motivating this approach. The ASED method proceeds from the view 
that chemical bonding is a process ensuing in physically distinct succes- 
sive stages, This is a theme that has some prospects of having theoretical 
validity. To illustrate this point, consider reviewing the simple case 
of bonding in the H 2 molecule-ion in the light of the ASED conceptual 
picture. Let the basis set for the variational energy calculation be a 
minimal basis set (MBS) consisting of a hydrogen atom orbital centered 
on each of the nuclei,^^^^ 

|A>=j~exp(- 

/ir 

jB> =^exp( - 

/it 


(7b) 
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The secular equation is 


«A.A - ^ 


^■A,B “ ^ ^A,B 


*^A,B ’’ ^^A,B ^A,A 


and the lowest-energy solution is 


( 16 ) 


(8) 


E “ XTTJ^ ^ ^A,B^ 

Starting from the usual electronic hamiltynlan for this caseOS), one 
obtains 



where I = -1/2 hartree « - 13,606 eV Is the first ionization potential 
of a hydrogen atom. One recognizes tho matrix element on the right- 
hand-side of Eq. (10) as the Ej^ energy of the ASED approach. Therefore, 
we write 


«a,a'-'"Er 01) 

In the case of the molecule-ion, Ep physically corresponds to the classical 
electrostatic energy for bringing a "rigid" and "non-exchanging" hydrogen 
atom up to a proton. Hence, we identify the ^ energy with the first 
process envisioned in the ASED conceptual picture of bonding. More- 
over, in simplifying the expression for Ep one obtains 




Er j P(r) - i) 


r^ dr 


(12) 
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where D » - J^gH and p(r) « 4 exp(- 2r). It 1s clear from this 

that ^0, and that, consistent with the ASED physical picture, E^ 

Is the energetically unfavorable factor for binding. 

The g matrix element Is tic f‘'plitude for the binding 
electron to make an exchange from basis \^> to basis [B>, or, equiva- 
lently, to "tunnel" from one atom to the other. Hence, we Identify 
qualitatively the g energy with the second process envisioned In the 
ASED conceptual picture of bonding. Starting from the usual electronic 
hamiltonlan for this case,^^®^ one obtains 

• - ‘ * ^''1 iTf%r " 

The evaluation of the matrix element on the right hand-side of Eq. (13) 

(15l 

shows that g g 0 for all except very small internuclear distances*' ' 
This Is consistent with the presumption of the ASED conceptual picture 
that the energy of the second process would be the energetically favor- 
able one for chemical bitiding. The identification of the role of ex- 
change matrix elements as the source of chemical binding, of course, has 
a long and niustrious history. Finally, these two energetically 
opposite factors, Eg and g combine directly in Eq. (9) to yield, as 
envisioned in the ASED conceptual picture, the approximate molecular 
total energy. There is in Eq. (9) a renormalization factor that does 
not affect these conclusions since it is a connion factor of both Eg and 
Ha g* The dependence on internuclear distance of the competing energies 
Eg and H^g, and the resultant energy, E, is illustrated in Fig. 1. The 
symbol A in the notation aE used in Fig. 1 indicates that the limiting 
value of E at large internuclear separations has been subtracted. We 
shall adopt this convention in the rest of this paper. 


Energy (hartree) 



L5 2.0 2.5 5.0 5.5 4.0 4.5 5.0 5.5 6.0 


Internuclear Oittonce(bohr) 

FIG. 1. The electrostatic energy E^, the exchange energy g, and the 
total energy change AE as a function of internuclear distance 
for the molecule in the X^lg state. 


■X 

I 
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The prospects of reconciling Anderson's ASED energy aiuatz with 
that of the MBS variational energy approximation are not good. The ASED 
energy for this case is 

E * E|^ + G {14 ) 

where 

«A.A ■= - ' 

’^A.B ' - 

These equations resemble those of Eq. (9) only in form. The irreducible 

difference lies in the use of the K factor in the ASED exchange matrix 

element, Eq. (17). This factor controls the contribution of the energy 

favorable for bonding to the total energy. Moreover, the total energy 

appears to be very sensitive to small variations in this factor. Figure 2 

shows the sensitivity of the ASED potential energy curve for Hg on the 

choice of the K factor. This figure also illustrates the large differences 

that can result between the MBS energy curve and the ASED energy curve. 

This is a numerical indication that the justification for Anderson's 

ASED energy formulae lies beyond the MBS energy formula. The exact 

solution may indeed be cast into the form of just a two-state problem by 

use of the partitioning technique. Such equations are the ultimate 

basis upon which the validity of the semi-empirical ASED energy formulae 

can be decided. The analysis of how semi -empirical approximations relate 

( 23 ) 

to the exact theoretical basis has been presented before by Freed.' ' 


11 


To conclude, Fig. 2 also shows, for purposes of comparison, the potential 
energy curve resulting from an exact numerical solution of the Hg molecule 
problem in the state. 

III. A FIRST^PRINCIPLES ASED METHOD 

We present below a different molecular-energy decomposition scheme 

than the traditional one.^^*^^^ One finds in this decomposition some 

(251 

common elements, however, to the work of Kitaura and Morokuma.' ' In 
the present work, we have attempted to adhere, as much as we find possible, 
to the ASED physical picture of the binding process to motivate the group- 
ing of the various molecular energy terms. 

The first task is to find the interaction energy associated with the 
process of building up the molecule by bringing together "rigid" and "non* 
exchanging" atoms from infinity to their assigned positions within the 
molecule. The interaction energy of this process is identified in the 
ASED energy aiucUz with Eq. (1), the Ej^ energy. Although this was found 
to be correct in the case of the Hg molecule reviewed above, it does not 
appear that Eq. (1) is, in general, the correct energy of this process. 

For instance, a first inspection of Eq. (1) shows that it contains, for 
each pair of atoms, an arbitrary choice as to which of the two atoms is to 
be chosen as the source of the electrostatic potential in this formula. 

It seems unphysical tiiat, in the absence of other factors, the correct 
interaction should contain an arbitrary choice of this type. 

It is possible to show that indeed Ej^ does not, in general, follow 
from the premises^^’^^^ from which it is proposed to be derived. 

Anderson starts the derivation of Eq. (1) by invoking the Hellmann- 
Feynman force-formula, 


Energy (10“' hortree) 



1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

Internudear Distance (bohr) 


FIG. 2. Comparison of various ASED potential energy curves for the Hg 

2 + 

molecule in the X Zg state with the minimum basis set variational 
energy^^^Vand the exact potential energy curve. 



12 


Va E 



% “ W 
li fe - fell ^ 



' fe^ 

lU - fell^ 



(18) 


where is the gradient operator with respect to the three components 
of the position vector and p{j(j) Is the exact molecular one-electron 
density. The exact density is then written in terms of the superposition 
of atomic densities and a complementary term, 


p(^) =E - y + 6p(j^) (19) 

In order to isolate the energy resulting solely from the superposition 
of the "rigid" atoms, Anderson advocates^^*^^*^®^ neglecting in Eq. 

(18) all terms not arising from the superposition of atomic densities 
in p(^). The result, which we denote E^, is 






^fe ' fe^ 

II Sb - ^aII ^ 



lU-Wl' 



(20) 


Note that, because each atomic density is of even parity under inversion 
through its own center. 



It - y 
II t ■ ssaIi^ 


dr = 0 


( 21 ) 


Using this result, one may write Eq. (20) in a more compact form 





(t - Ra) 

1 1 - SaII 


dr 

3 A/ 


( 22 ) 
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Where 


"B (j: ■ Sb> ■ h *(>: ■ - BB‘<t ■ Sb> (”> 

The function Ej^ Is obtained from the gradients of Eq. (02) by a line 
integral along a path taking the system from infinity to its assigned 
configuration in the molecule, 


Er -/S’a Er • “Ra 

This reduces, after some manipulation, to 



2 E 


A > B 



' 1 

IUa-%11 



124) 




II ^ / 


(25) 


The remaining integral over the x variable is a standard integral , 

The result is a complicated combination of terms the details of which 
become unimportant here. The point to note is that the summand in 
Eq. (25) is symmetric upon permutation of atom labels A and B. Thus, 
the problem with the arbitrary choice in the heretofore used formula 
for E|^, Eq. (1), is resolved. The corrected equation also shows that, 
if the atoms of a particular pair are identical, then the formula for 
that pair interaction reduces to just that obtained from Eq. (1). 

However, even the corrected expression for E^, Eq. (25), does 
not fully correspond to the energy of the physical process of building 
up the molecule from "rigid" and "non-exchanging" atoms. The source 


atiiii 


u 


of the problem Is that, when an approximate one-electron density, such 
as the superposition of atomic densities, Is used, then, the HeTlmann- 
Feynman formula, Eq. (18), only yields a part of the gradient of the 
energy. The other part of the gradient of the energy derives from 
the terms involving the partial derivatives of the wavefunctlon, There 
is in this case, no apparent physical justification for neglecting this 
other part. The total gradient, which we denote In this case Is 


Vet* 





’a "C 


f ■ y \ Sa> “i: ■'k'') 

where is given by Eq. (22), After some algebraic manipulation, one 
obtains 

'mjj - Ra) »b(i: - Rb) (37) 


where the density for example, is defined by Eq. (23). 

Proceeding just as before, one obtains the energy E^ from the gradient, 
® Integration, e.^j; Eq. (24). The result is 


'cZ 


_ i Y, E d);; d);;' 




II t “ x! II 


(28) 


This is just the energy expression expected from classical electro- 
statics. It Is tne work of assemblying "frozen*' charge distributions 
p^(r - R^) by the process of bringing each one successively from infinity 


~'V 
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to Its assigned position In the presence of the electric field of the 
"frozen" charge distributions already assembled, It therefore seems 
that I rather than Bp, best satisfies the requirements of the correct 
energy expression for the first stage of bonding of the ASED approach, 
The next step Is to Incorporate the E^ energy Into the molecular 
energy. The expression for the molecular total energy 


E * E g^ + 2-» 2-r 


hh 


lUft-gull 

The exact electronic energy may be written as the expectation value 


(29) 


' ‘''^1 (30) 

where, as usual, the exact many-electron molecular hamiltonian, H, 
is the sum, over all electrons, of the kinetic-energy, nuclear- 
attraction, and electron-repulsion operators^^®*^^^ 


fi.ihCi) E-~-! — 

i 1 > J lUi - rjii 




(31a) 


B IU-SaI1 

The exact many-electron wavefunction, T, is one of the eigenfunctions 
of the hamiltonian operator, H. The particular electronic state to 
which refers to is dictated by the choice of the eigenfunction T. 
Consider constructing the following combination of terms, which we 
label D, 


(31b) 
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D m 


iT 


A B 


( 2 

ik-W 



~ 8b> 

IIc-SaII 




tff "B ' 8b) (32J 

JJ lU - );* II 

Adding and subtracting the term D from the right-side of Iq, (29) yields 
the molecular total energy in the form 


^ • ^o£ + ^ee “ ° ^33) 

where is defined by Eq. (28), As discussed above, E^ is the term 
that physically is identified with the first-stage of bonding in the 
ASED approach. 

In order to show how the D term combines with the term, we 
construct a convenient energy decomposition scheme. We partition the 

a/qN 

many-electron hamiltonian into a zero-order hamiltonian» H' % and a 

/s/l\ 

perturbation, \ 


h.h«>Wh(’) 

where the terms of the partitioned hamiltonian are defined as 

H<'».L(h(i) + 0(i)) 
and 


H<1) -L E 

i > j 


1 

Ui ■ iCjll 


-E?(i) 


(34) 


(36) 


(36) 
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This type of partition Is well-known trom perturbation- theory approaches 

to the electron-correlation problem/^^"^^^ In the electron-correlation 

problems* v Is usually chosen as some form of the Hartree-Fock potential. 
For a closed-shell system, for example, this 

where the Coulomb and exchange operators are defined respectively 

The summation ranges over all occupied Hartree-Fock molecular orbitals, 
t|»|^ ()r). In the present case, consistent with the physical picture of 
the first-stage of bonding of the ASED approach, we define v Instead as 
a superposition of atomic potentials, each centered on its own atomic 
nucleus, 

V = E vf (38) 

A " 

There Is no unique choice of the *, /Omic one-electron potentials 
v^*. Physically, we wish these potentials to enable an electron of 
the molecule, when In the vicinity of one of the atoms, to interact 
with the electrons of this atom as though it were part of the electronic 
structure of this atom.’^®^ In the Hartree-Fock approximation, the 
potential meeting this physical criterion (in a closed-shell system) 
is given by Eq, (37) but as constructed from atomic-orbitals rather 
than molecular-orbitals. In more accurate approximations, a similar 
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atomic potential exists but its form may be much more complex. In 
general, an approximation for an n^-^electron atomic system which 
retains an independent-particle interpretation will have at most n^^ 
different atomic orbitals. Each of these different atomic orbitals 
can be interpreted as the state of an electron moving in the field 
due to the other electrons. That is, each atomic orbital , <l>|^ is 
the solution of an equation 



“•A 


lU - 8 a1I 




'•'k.A ” 'k,A '•'k.A 


(39) 


where 0|^ ^ contains the interaction effects due to the other 
electrons of the atom. Each occupied orbital is associated with its 

A 

own characteristic (n-1) -electron potential To accommodate 

the specificity of each potential to its atomic orbital, we write, in 
general. 


''a*'’ ^'^k,A l'*’k,A><''’k,Al 


(40) 


where we implicitly assume in writing the equations in this form that 
the atomic orbitals, i|)|^ form an orthonormal set of functions. In 
the case that these orbitals are non-orthogonal , a somewhat more com- 
plicated expression applies. Note that the summation in Eq. (40) is, 
in principle, not restricted to just the n^ occupied atomic orbitals. 

In practice, however, one may obtain a manifold of physically meaning- 
ful unoccupied levels from the potentials for the energetically-highest 
occupied orbitals. For example, one may in such a case write 
Eq. (40) as 
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The partition of the electronic haiiHltonlan Into a ^ero^order 
hamlltonlan, and a perturbation, are the basic Ingredients de- 
fining a perturbation expansion of the electronic energy, Eq. (30), os follows 


’qZ 


E^fl^ + + E^P + E^P + . 

^cZ '^cZ 

:( 0 ) 


(42) 


The zero-order energy, E^^ , Is the eigenvalue of the zero-order 
hatnlltonian An eigenfunction of Is a product of N one- 


electron functions 


'*’^^Cl,)^2, • • •’ '*'2^*C2^ * ‘ ' 

where the one-electron functions, >|)|^(j;;), are eigenfunctions of the 
equation 

R II K \ II 

N is the number of electrons In the molecule. The actual distinct 
number of occupied orbitals appearing In the basic product wave- 
function, Eq. (43), Is determined by the usual building-up principle 
of occupying the energy levels, C|^, in order of Increasing energy and 
in accordance with the Pauli-principle. The zero-order energy is thus 
the sum of the one-electron molecular orbital energies 


^ ^%cc^k 

where (n^--)r, is the occupation number of the k'th molecular orbital 

OCC K 

level, and the sunination ranges only over the occupied levels. Al- 

*(ol . 

though ^ is an eigenfunction of H' f, any permutation of the 
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electron-coordinates also produces a degenerate eigenfunction. The 

A / Q \ 

physically most complete eigenfunction of H' ' is therefore, an anti- 
symmetric combination of these products which is an eigenfunction of 
spin,^^^^ i.e. 

= A 1^4. (L C. (46) 

S N 

where A is the antisymmetrizer, and 0 ^.’ is a linear combination of 

N-electron spin-function products which yield 

spin 

= S(S + (47) 


we indicate a sum over spin-eigenfunctions 0 ^*^ in Eq. (46) since, 
as is well known, there are generally several linearly independent 
spin eigenfunctions of spin S arising from the coupling of N spin 
1/2 particles. 

Expanding the molecular-orbitals, in terms of a basis set 
of atom-centered functions, x . a 


= \,A<C - «A^.A;k 


(48) 


leads from Eq. (44) to a matrix-equation such as Eq. (2). The 
hamiltonian matrix-elements are defined as 


P,A jv ,B 



(49) 


The evaluation of the matrix elements of the usual one-electron 
operators poses no particular problem. The matrix elements of the 

Aa ^ 

atomic one-electron potential, v^ , may be evaluated conveniently by 
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preparing, from the atomic orbital basis, a matrix- representation of 
the short-range part of the potential. For example, define an effective 
charge as 

- Ua) = ("a - d [i - ®xp(- “1^ - &l^)] 


where n^ is the number of electrons on the A'th atom. The parameter 
a is defined so that the onset of the asymptotic behavior of 

Zf^(r) 'v- (n^ - 1) (51a) 


A 

coincides with the onset of the asymptotic behavior of 



(51b) 


The difference potential, 


A ■ Ra> 


(52) 


is then defined as the short-range part of the ^ potential. This 
short-range potential may be represented in a basis centered only on 
the origin of the potential. We write this representation as 


n fA!LJA) .y y 




(53) 


where 






dr 




(54) 


The convergence of this representation is expected to be rapid given 

A 

the generally slowly varying character that is typical of Uj^ ^ potentials 
in this inner range. Therefore, we expect that restriction of the 
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representation to a finite basis yields a reliable approximation of 
the short-range properties of the potential. The net expression for 
the hamiltonian matrix elements is 




<x 


y,A 





( ” tr ) 

lnF¥ii 



(55) 


If the basis functions x aare chosen as the atomic orbitals, then 
the expression for the hamiltonian matrix elements takes on yet another 
simple form. We find 


H 


VI 


,A;v,A 





(56) 


and 
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y,A;v,B 


_ 1 


i (e 


p,A 




y ,A;\j,B 
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y,A 
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IIV - Sa« " 



C/B 

The one-electron energies, e ., are approximations the true 
ionization potentials of the atoms.' ' Hence, through Eqs. (56) 
and (57), the hamiltonian matrix has a superficial resemblance to 
the corresponding matrix in the ASED approach as defined by Eqs. 
(3a) and (3b). If the one-electron atomic potential, Eq. (40), can 
be reduced to one common potential for all atomic orbitals, then 
one obtains, even for an arbitrary basis, the simple result 



(■ 


? + 


) 
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+ < 




(58) 


Any of these formulae are very convenient, from a practical point 
of view, because their calculation, we note, requires at most the 
evaluation of one -electron- type integrals. 

Thus far we have established the basis for a convenient decom- 
position of the exact molecular total energy into the form 

We find that the first two terms, have a resemblance to 

the semi -empirical ASED energy formulae, Eq. (6). We consider next 
the terms in Eq. (59) that are required to improve the approximation. 
In particular, we consider the nature of the first-order energy and 
its combination with the D term. The first-order energy is defined 
as 


eO) .<,(o)|j)0)|,(o)>^<,(o)|^{o)> 


(60) 


This energy may be written generally as (22) 


A# ‘ 




where is the diagonal component of the two-electron 

density matrix after integration over the spin variables(^^^, and 
is the one-electron density matrix, also after integrati 


on 


r 


E I 
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(22) 


over the spin variables.'^ ' Note that we shall adopt the convention to 
denote the diagonal component of the one-electron density-matrix, namely 
as simply p(r). This is also, in fact, just the electron- 
density function we have been using in the foregoing discussions. To 
proceed, we find it convenient to designate a part of the two-electron 
integral in Eq. (61) as a Coulomb energy, and the remainder, as an ex- 
change energy. For this purpose, we use the device of defining formally 

(39l 

a two-electron correlation function as' ^ 

In terms of this definition, the first-order energy assumes the con- 
ceptually convenient form 


V 2JJ\\ 


p(r)p(r. ) 


n 

r - r 


drdr 
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xc/t 




where 
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dr dr' 
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The D term, Eq. (32), may also be written in the convenient form 

at/ 
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M- 


<>(r) - 


sup 


r - R„ 
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(63) 


(64) 
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where Is defined as the superposition of atomic densities, 


-s Pf(. - s,) 

'U n 

Combining Eqs. (63) and (65) to obtain the expression for (E^^^-D), 
we find 

Z, 


(66) 


A + (;at 






■?J '■ 

l|r - ir'||‘’A*X"SA*‘’r**:'‘SA^ ‘‘t “'t' 

Haying completed the definition of (E^P - D), we seek to 


( 0 ) 


(67) 


combine the energy terms of the zero-order energy, > with the terms 
(E«P - D). First we note that E^P, Eq. (45), may also be written 


as 


■qZ 






x: 
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Combining this formula for with D), Eq. (67), results in 

the cancellation of some terms. The net resulting formula Is 


:(0) , ,(1) 


*ee 


u 




^ ^ 4c/i 


( 69 ) 


where 


A/ "'V 

and 

«p'^’(c)«p‘”(C) <>K <<«' (71) 

and 



Equations (70) and (71) give the effect on the energy of the change 
in the one-electron density matrix relative to the superposition 
of atomic densities. The relaxation in the density matrix is de- 
fined as 


6p(l^(r,^') = P^^^(r.r') - P^^P(r,j(;;' ) 


. 0 ) 


(73) 


Note that 




dr = 0 


(74) 


:(0+l) 


The term ■' represents the linear effect on the kinetic, electron- 

nuclear attraction, and Coulomb electron-electron interaction energies 


r 

r 


[ 


11 

due to relaxation of the electron density. The term the 

non-linear relaxation effect. Both and 

the same as E^, approach zero as the atoms separate to the neutral- 
atoms limit. The E^ energy is constant with changes in molecular 
geometry. It is the sum of the kinetic, electron-nuclear attraction, 
and Coulomb electron-electron repulsion energies of the isolated atoms. 
The formulation of the sum (E^^ + E^^) remains incomplete with- 
out a more precise definition of the exchange energy, For this 

purpose it is necessary to specify in greater detail the structure of 
the zero-order wavefunction Eq, (46), Let the orbital product 
<I> in the zero-order wavefunction be composed of n doubly occupied 
orbitals and m singly occupied orbitals, 


* ' ’^M**,'^n Vl V2 **• Vm 

where 2n+m = N. The one-electron density matrix resulting from this 
wavefunction, Eqs, (46) and (75), is 

n m 


(75) 




k=l 


(76) 


k=1 


The exchange energy associated with this wavefunction is 


(40) 
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1^1 ^ 1^1 ^1 2 ^k£ *^k+n,^+n (77) 


where 


hi ^i 


= 1 




k=-e 


(78a) 


(78b) 
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The in Eq. (70) are just the exchange Integrals defined, as usual, 
by 


.rrhWK> 

Kici ’Jj nipyii — 


<<c' 


(79) 


All the properties pertaining to the total spin of the system, S, 
and to the relative weight of each alternative linearly independent 
spin coupling of the m unpaired orbitals are contained in the gj^^ 
coefficient. The number of these linearly-independent spin- 
couplings is denoted by n . The value of n is determined by the 
desired total spin and the number, m, of unpaired orbitals. It is 
given by^^^^ 

n « (IS , t. . (80) 

" (^m + S + 1)! {-im - S)! 


The properties of the gj^^ coefficients have been discussed extensively 
by Goddard, Ladner, and Bobrowicz, The c^. in the expression 

for gj^£ give the relative contribution in the total wavef unction of 
each alternative spin-coupling. They are obtained from the eigenvectors 
of the secular equation 


l[( 


k> I 


^•^k+n, £+n)“ \ ^ij] ‘^j;k “ 


0 i,k«l,...,n- (81a) 


The eigenvalue A|^ is the energy difference 


* 

a 


(81b) 


★ 

The energy is independent of the details due to spin and spin- 
coupling. It is only dependent on the number of doubly and singly 
occupied orbitals in orbital product, Eq. (75). It is given by^^®^ 


# n A ni A n n 

m n I m m 

■*■ Jij ^^'^k+n,^ " ^ $\ P\ '^^+n,^+n ( 82 ) 

where the one-electron operator n Is defined 1n Eq. (31b) i and the 
Coulomb Integrals are defined, as usual, by 


k,.t 


rr fklt) 4(«') j,. 

’Jj Ik - i'TT 


(83) 


The average energy of the n. possible spin-coupling states follows 

V- 

from Eqs. (81) and (82). It Is 


(eUS’ + E<1>) 


(Ll 'avg 


ht* 


m m 




where 


V) = i <v )/"s 


(84b) 
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All the information about total spin and spin coupling affecting the 
total energy resides In the uJj"’(T,,p matrix elements. These are the 
elements of orthogonal matrices that yield the Irreducible repre- 
sentations of the symmetric group, based on standard Young 
shapes and tableux, The matrix elements used here belong 
to the matrix representation of the elementary transpositions t|^£. 
x^’'’^(Tk£) of Eq. (84) is, apart from the n^ divisor, just the 
character of the (S,m) Irreducible representation for the element 


The exchange energy can also, at least formally, be decomposed 
into a term constructed from the limiting orbitals to which the 
molecular orbitals tend as the atoms separate, and terms that 
give the relaxation effects relative to the first term. 


■ E° + 6e1I1 + 


•xc/i ’ ‘^KCh ^ ” ‘^xcJt 

Let the set of molecular orbitals become the set of limit orbitals 
{i|i^ las the atoms separate. Then given by Eq, (77) but 

replacing the orbitals with the limit orbitals to construct the 
various terms. The relaxation energy in turn is 


(85) 


i / ( • ll ' I (8«) 


where 


«p^’>(r,r') - (ho„)^(fk(j)tKtt’) - 4'5(j;)tk(]C'» (8^) 

and 

<Vc>k " 2 Kk<n (88a) 

“ 1 n+l<l«n+m (88b) 

The i<£ are just the usual exchange operators, Eq. (37c). The super- 
script is used here to indicate that the orbitals used to construct 
are the limit orbitals The non-linear relaxation effect on 

the energy is 
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Combining the expressions for the various relaxation energies, 
Eqs. (70), (71), (87) and (90), one obtains the following expression 
for the ASEO energy through first-order 


c(0+1) * p8t . p . pO + 

^ASED * -U hdi ^^ASED ^ ^ASED 

Moreover, the energy term linear In the density-matrix relaxation 
can be written as the expectation value of Fock-type operators, 
Combining Eqs. (70) and (86) one obtains 




ksi 
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AQ AQ 

where the Fock operators F and F|^ are 
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(90) 


(91) 


(92a) 

(92b) 


When the limit orbitals, are each just the atomic orbitals, 

as in the case of all heteronuclear molecules, the Coulomb and 
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exchange potentials may be grouped according to the atomic nucleus 
of the atomic orbitals from which the potentials are derived. In 
those cases, the Fock operators consist, apart from the kinetic energy 
operator, of superpositions of atomic effective potentials 




(93a) 


F? = - 59^ + E V? 


k 


A 

where the atomic effective potentials are defined as 

Z, 


(93b) 




1 r>o 


A 




“ - rtn ^ £c ^ £0 ^ 


C 0 

The and denvlifi respectively the sets of integer labels of zue 


doubly and singly occupied atomic orbitals of nuclear center A. 

The net non-linear effect on the energy due to the relaxation in the 
density-matrix is obtained by combining Eqs. (71) and (89). One 
obtains 




+ E E ( 29 t,p + 1)«P*'*(r,r')Spr)(r,r')ldl- dr' 
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( 95 ) 
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The correlation energy has been defined previously, Eq* (42), as 


pcorr » 55 r(2) ^ p(3) ^ 


( 96 ) 


This perturbation expansion of the electronic energy follows along 
conventional lines^^^”^^^ from the definition of a partition of 
the electronic hamiltonian, Eqs. (34) -(36), and Eq. (38), The 
second-order Rayleigh-Schrodinger perturbation correction to the energy, 
for example, is 


e'c* = - E E Kv 

^ k a 


n(0) I iTi ,„avi 2 




The and are the configuration functions generated, re- 
spectively, by replacing one or two of the occupied orbitals 
(denoted by k,£,...) in the orbital product of the zero-order wave- 
function with unoccupied orbitals (denoted by a,b,..,). The 
second-order energy, E^% has non-vanishing contributions from 
single-excitation configuration functions because the molecular 
orbitals used here do not satisfy the Brillouin condition. 

This perturbation formulation is not directly applicable when 
an incorrect admixture of covalent and ionic terms in the wavefunction 
causes extraneous long-range terms to apijaar in the first-order energy. 
The terms in the correlation energy that cancel this extraneous long- 
range behavior derive, as is well known, from the configuration- 
functions involving excitations to those unoccupied molecular orbitals 
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that become degenerate with the occupied orbitals as the atoms 
separate. In order to continuously treat the correlation correction 
to the exchange energy, a prior transformation of the hamiltonian 
matrix has to be made that properly resolves the increasing de- 
generacy in the domain of large internuclear separations. One method 
applied to a variety of problems of this type is the Van Vleck trans- 
formation. Consider the hamiltonian submatrix that derives from 
the subspace consisting of and all configuration functions gen- 
erated by excitations to unoccupied orbitals that, as the atoms 
separate, become degenerate with the occupied orbitals. A straight- 
forward application of the Van Vleck transformation modifies the 
coupling of the sub-matrix to the remaining matrix so that the modified 
off-diagonal coupling terms are reduced to be of second-order or higher 
in the perturbation. As is well known, the eigenvalues of the actual 
transformed sub-matrix are then automatically accurate through third- 
order in the perturbation expansion, since the modified off-diagonal 
coupling terms can only contribute to these energies in fourth-order 
or higher. 

The correlation energy can also be obtained by the more tradi- 

(44) 

tional variational configuration interaction (Cl) approach.' The 
standard Cl variational energy is 

? ‘^i ^ij 


(98a) 
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where 


HlJ {98b) 

and ' 1 '^ stands for the 1'th configuration-function, and 'I'l is defined 
as the zero-order wavefunction, Each diagonal haniiltonian 

matrix element can be decomposed as shown in Eq. (91). Note that 
the latter equation is just a prototype for all diagonal matrix- 
elements since 

Each element, in general, has its own particular separated-atoms 
limit characteristics relative to \ifhich the relaxation terms are 
defined. Extending these arguments to the off-diagonal haniiltonian 
matrix elements, we write formally 

Hii = ( Eot)i1 ^ 

"lj “ "ij * noob) 

where is the off-diagonal element constructed from the limiting 
orbitals of the configuration functions 4'^ and 4'^, and 6H . . and 
are the relaxation effects relative to h9.. 

I J I J 
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IV. DISCUSSION 

The molecular orbitals energies, obtained by solving the 
ASED one-electron equation, Eq, (44), become the atomic orbital 
energies in the limit as the atoms of the molecule are infinitely 
separated. The reason is that, in this limit, the ASED orbital 
equation describes a single electron in the field of the isolated 
atomic potentials. The resulting spectrum of one-electron states 
in this limit is, therefore, just the superposition of the one- 
electron states of each of the atomic potentials. This property 
enables one to follow in a continuous fashion the evolution of the 
ASED molecular orbital energy levels from their parent atomic levels. 
Moreover, in this limit the molecular orbitals, \|)|^, become either 
the atomic orbitals themselves, or, otherwise, just spatial -symmetry- 
dictated combinations of these atomic orbitals. 

The ASED zero-order energy, Eq. (45), is determined by occupying 
the molecular orbital levels resulting from Eq. (44) in order of in- 
creasing energy, e|^, according to the usual building-up principles. 

This choice of the orbital occupation may not always lead to the desired 
diabatic electronic state at all internuclear separations. A different 
orbital occupation may be characteristic of the adiabatic state in these 
other domains of the energy surface. In these cases it clearly becomes 
necessary to go beyond the ASED energy through first-order, Eq. (91), 
to a configuration interaction (Cl) calculation of the energy. The con- 
figuration functions are the various zero-order wavefunctions, Eq. (46), 
corresponding to the different choices of orbital occupation. Despite 
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the increased complexity, the ASED energy decomposition scheme can still 
be used even in this case to simplify the Cl energy calculation. 

A. The LiH molecule 

To illustrate the above, we consider formal aspects of the ASED 
application to the LiH molecule in the state. The molecular orbitals 
obtained by solving Eq, (44), {I 0 , 2o, 3o, 40 ,}..., approach the limits 


10 - 


(101a) 

20 — 

— Uh 

(101b) 

30 — 

— 2sli 

(101c) 

40 - 

— 2p»|_, 

(lOld) 


as the Li and H atoms separate. The application of the building-up 
principle leads to the following ASED zero-order wavefunetion 




= A (lo^ 20 ^ a3 a3) 


( 102 ) 


As the atoms separate, this wavefunetion approaches the Li^ (^S) 
+ H"(^S) limit. The one-electron density matrix deriving from 

^^( 0 ) ^5 


= ZiZo.Za) (t03a) 

At large internuclear distances, approaches the superposition 

of atomic-ions density matrices 




(103b) 
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We adopt here the convention of using a dot between each pair of 
orbitals in the expansion of the density matrix to indicate that they 

depend separately on the spatial coordinates jjj; and In contrast 
the adiabatic electronic state approaches the neutral atom 
limit Li^^S) + H(^S) at large internuclear separations. The ASED 
zero-order wavefunction describing this limit is 

4^^ - A ( la^ 2a3o aB (104) 

With obviously different orbital occupations than obtained by appli- 
cation of the building-up principle. The one-electron density matrix 
deriving from 4^^ is 

= 2(lo.lo)+ 2o.2a + 3a. 3a (105a) 

At large internuclear distances, r') approaches the super- 

position of neutral atom density matrices 

The Cl wavefunction is a linear combination of the two alternative 
zero-order wavef unctions, Eqs. (102) and (104). The diagonal 
elements of the hamiltonian matrix in this configuration function 
space are the same as the energies through first-order for each zero- 
order wavefunction separately. The ASED decomposition for each 
diagonal hamiltonian matrix element is obtained in Eq. (lOOa). 
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Whereas the ASED decomposition for the ^22 element Is based on the 
superposition of neutral atoms In the first stage of bonding, the ASED 
decomposition of the H-ji element Is based Instead on the superposition 
of atomic Ions. Thus, at large Internuclear separations, D. 
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while (E^)g 2 falls off exponentially to zero. The detailed 
formulae for (^Jc,/t^‘1 special cases, respectively, 

of Eqs. (28), (72) and (77). We divide (eJ^^)^^ Into those terms that are 
atomic exchange energies, and hence remain constant, and those terms 
that vary with Internuclear separation. We write 




In the present example, these terms are 
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*^1Sli,2Sl^. 


■2Sl, ,2S|_j 
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2 ^ls,^ls„ 

(108c) 


"KiSh,2Sl,. 

(108d) 






The general formulae giving the effects on the energy of relax- 
ations in the one-electron density matrix are Eqs. {92), (93) and (96)* 
For purposes of illustration, we next give just the linear relax- 
ation terms for this case. First, relative to the superposition 
of atomic-ion density matrices, one obtains 

= 2 <lo - Is, jFllo + Is, .> + 2<2o - lSnlF^l2a + TSn> (109a) 


where 


A P A / A A \ / ^ ^ 


Second, relative to the superposition of neutral -atom density matrices, 
one obtains 

SHjj = 2<lo - lS|_, + lS|_^> + <2o - lS(,|F^^|2o + ls^> 


+ <3c - ^ 


where 
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Both the and terms are seen to fall off exponentially 

to zero at large Internuclear separations. The long-range Inter- 
actions, If present, are carried by the terms. The net 

limiting values of at the separated-atoms limit Is given by com- 
bining with one obtains 


and 

(e'J ^ ^>22 = 'Ll ^ Eh 
+ •• 

E^^ 4 . and E|^, are the Li and H Ion energies In the Hartree-Fock 
approximation, 

Eli+ « 2<lSL^ihL^iiSL^> (^'^ISl^.ISl^ " '^''^Lr^^l) 

Eh- = 2<ls^|h^lls^> + 

E|_^ and E^ are the. analogous energies for the Li and H atoms, 


(nia) 


(n2a) 

(112b) 


Eli = K2sLil\li2^L1> * (eJisl^.Es,, ' 

E[^ - I b|^ i 1 (113b) 

A A 

In these equations and h|^ are the one-electron operators, Eq. 

31b), for the Li and H atom cases respectively. 

The off-diagonal hamiltonian matrix element Is folded into a 
form analogous to that obtained for the diagonal elements as indicated 
in Eq. (100b). 
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The first term of Eqt (100b) in 1 
where 

fc . h t 

In relaxation terms are 
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" 2 (la.lo - ISl^.ISl^) 
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B. The Hg molecule 

A prototype case Illustrating a quite different challenge for 
the ASED energy decomposition Is found In the Hg molecule In the 

Xl Sg state. The ASED orbital equation In this case Is the same 

y 

as that of the Hg molecule-ion problem. The molecular orbitals 
obtained from solving the equation, {lOg, 10 ^,...}, approach the 
limits 
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as the H atoms separate. The two hydrogen atom orbitals are labeled 
1s^ and Isg, The ASED orbital energies, however, both approach the 


ls» as 

the atoms separate 
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The application of the building-up principle leads to the following ASED 
zero-order wavefunction 

4'^°^ « A (lo ^ a$) (119) 

The zero-order energy is 

fr\\ 

(120) 


:( 0 ) 
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= 2 


It follows that the zero-order energy approaches the correct limiting 
energy as the atoms separate. 
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However, at large Internuclear separations, the zero-order wavefunction, 
Eqn. (119a), becomes 



[( 




j o6 + A + ISgjaS 


The first component correctly describes the + H(%) limit. The 
second component of this wavefunction, however, describes an ionic 
limit. This admixture of the covalent and ionic wavefunctions even in 
the separated-atoms limit is the cause of extraneous long-range terms 
in the first-order ASED energy to be discussed later. 

The zero-order ASED energy is defined as 


p(0) Bp + p(^) 

^ASED V 


( 121 ) 
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The two components of E^g^Q superficially resemble the E^^ and 
components of the energy cmatz in Anderson’s semi -empirical ASED 
method, Eq. (6), Figure 3 shows that E^, except for a small attractive 
region, is mostly a repulsive energy curve. Figure 3 also shows that 
aE^^, in contrast to E^, is a strongly attractive energy curve. The 
nature of these curves is consistent with Anderson's qualitative argu- 
ments. In these calculations we have used the approximation of re- 
stricting the orbitals to a minimum basis set In order to simplify the 
calculations. The net energy curve, AE^g|g, is also shown in Fig. 3 
together, for reference purposes, with the exact potential energy curve. 
This comparison shows that the zero-order ASED energy is quantitatively 
in gross error. The binding energy is overestimated by an order of 
magnitude. 




Internucleor Distance (bohr) 


FIG. 3. The zero-order eleitronic energy the classical electrostatic 

energy E^, and the net ASED zero-order potential energy, 
as a function of internuclear distance for the tig molecule in the 
X^Eg state. The exact potential energy curve^^®^ is also shown 
for purposes of comparison. 
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The first-order ASED energy is defined as 


4seo’ f 


The linear relaxation energy Is divided Into contributions from the kinetic, 
electron-nucleus, and electron-electron repulsion energy 

The various relaxation energies are defined as 
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where the relaxa-^ion in the one-electron density matrix Is 

= 2(lo. lo.)- 


«p' 
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The non-linear effect on the energy due to the relaxation in the 
electron density is 
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Note that we have dispensed with the division of the electron- 
electron repulsion energy into a coulomb and an exchange contribution. 
As is well known, in this case, the first-order energy contains terms 
that cause this energy to approach the wrong limit as the atoms sepa- 


rate. 


( 47 ) 


The partition of the electron-electron repulsion energy only 
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serves to obscure the problem in this case. The terms in the ASED 
energy, Eq, (122b), that cause this incorrect approach to the limit 
are isolated in the term. This term is 




p(0+l) 

"A '“A, 

with the property that, at large internuclear distances, D, it be- 
haves as 


( 127 ) 


r(0+l) 




^ ('^’=A.''=A ■ 0 


(128) 


A ^ 

The first term in Eq. (122), contains the correct separated - 


atoms limit energy, namely 


^el '' ^®ls 


(129) 


The various relaxation energies are shown in Fig. 4 as a function 
of Internuclear distance. This figure shows that the relaxation 
effects in the kinetic energy and the electron-nuclear attraction 


:(0+l) 


,( 0 + 1 ) 


energy, ' and ^ respectively, are by far the dominant 


contributions to the total relaxation energy. The linear relaxation 


effect in the electron-electron repulsion energy, is shown in 


Fig. 4 to be a comparatively small and slowly varying. The non-linear 
relaxation effect, is shown in Fig. 4 to be negligibly small. 


ae 


:(0+l) 


The total linear relaxation energy, 6E^ V, and the classical 


electrostatic energy, are shown in Fig. 5 as a function of inter- 


nuclear distance. As already seen on Fig. 3, E^^ does have a small 
attractive region. The linear relaxation energy, however, is the 
dominant attractive energy contribution. At small internuclear distances. 


I >if Mi liffil I jjiirtlililll'jl il I i f III 1 1 liiilii li n H liT . - 


Internuclear Distance (bohr) 


FIG. 4. The relaxation in the various components of the electronic energy 
as a function of internuclear distance, for the Hg molecule in the 
state. 





Energy (I0"‘ hartrae) 



FIG. 5. Decomposition of the partial ASED energy through first-order, 

^^ASED^ - into the classical electrostatic energy, E^ , 

and the relaxation energy, 6 e|J''’^^, for the Hg molecule in the 
state. 


dt al- I 
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the repulsive energy Is derived dominantly from The term 

is not by itself a meaningful term. We show later that there are terms 
in the correlation energy that precisely match this term, and that, when 
combined with yield a net meaningful interaction. For this reason, 

we only show a part of the first-order ASED energy in Fig. 5, namely 
(^^ASED^ " This partial energy is in fact just the sum of E^^ 


and 6E 


( 0 + 1 ) 


, The A in the notation AE^g^J^ indicates that the energy 
of the separated atom, E®^, has been subtracted from E^sED^' ® 

shows an alternative interesting decomposition of the partial energy 
(^*^ASED^ “ figure we show that, when and 

are combined with E^|^, a purely repulsive effective electrostatic 
type of energy is obtained. This decomposition focuses attention on 
the relaxation in the kinetic energy, as source of the attractive 

energy leading to chemical binding.^^^’^*^’^®^ Figure 7 shows a comparison 
of the partial first-order ASED energy curve with the exact energy curve^^^^ 
and the Hartree-Fock energy curve. We note that there is a consider- 
able improvement in the comparison with the exact result in going from 
the zero-order curves. Fig. 3, to the present one. The binding energy, 
for example, is now within 5« of the exact result instead of being off by 
an order of magnitude. Detailed properties of the partial first-order 
energy alone, however, are not in good agreement. The equilibrium bond 
length is overestimated by about 0.5 bohr, and the range of significant 
binding extends to internuclear distances that are much too large. The 
Hartree-Fock energy, E||p, is shown in Fig. 7 for comparison purposes in 


the form of a Nationalized" potential energy curve. 


(49) 


namely 


V. 
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Internuclear Distance (bohr) 


FIG. 6. Decomposition of the partial ASED energy through first-order, 

into an effective electrostatic energy, 

^cZ ^ relaxation in the kinetic energy 

for the Hg nrolecule in the state. 
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FIG. 7. Comparison of the partial ASED energy through first-order, 
^^ASED^ - with the Hartree-Fock^^^^ and the exact^^ 

potential energy curves, for the H2 molecule in the 
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- E^ (130) 

The first-order ASED potential energy curve is equal to this 
“rationalized" Hartree-Fock potential energy curve, i.e. 

^^ASED^ 

except for negligible second-order effects. It appears from Fig, 7 
that either aE^^^ or (aE^j^q - are equally unsatisfactory by 

themsel ves . 

As alluded to previously, the term has a direct counter- 

part among the terms in the correlation energy. Let us consider the 
configuration interaction (CI) approach to the calculation of the 
correlation energy. The configuration functions^^®^ are 

4'^ = A(lOg^aB) (132) 

and 

H'2 = A(la^aB) (133) 

where is just the zero-order wavef unction, Eq. (119a), and the 
la^ orbital in 4*2 is the unoccupied solution of Eq. (44). The total CI 
wavef unction is 

4' = 4'^ + C 2 4*2 (134) 

and the corresponding variational energy is 



( 135 ) 
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The hamlTtonian matrix elements can each be partitioned as we have 
already indicated in Eq. (100). The ASED molecular total enerfly, 
Eq. (29), is then obtained as 


- pSt X p X 

^ASED “ ^c£ 


^^1 ^^11 


:(0+l) + 


6H 


12 



+ (136) 

p 

where 5 E contains all effects that are second-order in the relaxation 
of the appropriate one-electron density matrices. The relaxation 
matrix elements, are divided again into contributions from the 
kinetic, electron-nuclear attraction, and electron-electron repulsion 
energy, 





(137) 



(138) 


These three types of terms are defined as 



•‘1C 039a) 
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(’39c) 

The relaxation In the density matrices of each Cl matrix element is defined 


as follows: 

<Spii^ ():<') * 2<l0g.l0g)- (T40a) 

* 2(l0y.lo^)- (HOb) 

'^‘^iP^Ju’^' ^ (140c) 

and the superposition density matrices in this case are 

o?r(w') = 9'“^ (*:•):') (’-”») 

9?r («’) = ISa-’^A- 


We note that, the same as for e[P), Eqs. (90) -(94), the 

relaxation matrix-elements, Eqs. (138) and (139), are clearly integrals 
of just atomic operators acting on relaxation density matrices. 

The molecular total energy, Eq. (137), can thus be written as 


-Cl 

*ASED 




where the correlated linear relaxation energy is 


5E^» = E?^. + 6E^^ + 6E^^ 
q£ fun m cc 


(143) 


The individual types of energy contributing to the linear relaxation energy 
are 

rr /' / ^ r\ /i\ 

(144a) 




^ *v 
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en 


/(h" 


6a II lU - 1 




6E 


Cl 

ce 


ifh^ 


+ .sup 


^12 


II i; 
<6) (' 




2 *Pi2'l6') *'z) 


. 0 ), 


K ''6' 


(144c) 


where the relaxation in the Cl one-electron density-matrix is defined as 


spci’lW'l ■ '■ci'lK'l ■ p^“’’IW> 


. 0 ) 


sup I 


(145a) 


and 

Pci’<W'> = ‘=1^(2 ’"g-’”g)* ‘=2^ 

A1 so , 

p(corr) , 2 C^ Cg 

The Cl coefficients C^ and Cg are the components of the eigenvectors 
of the secular equation that follows from Eq. (137). Neglecting non- 
linear relaxation effects, these Cl coefficients are given by 

C, = (f<°+’> + 6H,2y^ (6H„ - t.f + (147a) 


Cj ■= - (sH,, - (SH,, - 4)“^ + (147b) 

and the Cl eigenvalues are given by 

A = ^ («H,, + SHjj) ± ^ I (SH,, - + (F<°^'l + (>'•8) 


"".■'■vir-r; 


As the separated-atoms limit is approached, the Cl coefficients 
of the lowest energy root approach the limits 

C, (149a) 

C :1 (149b) 

and, hence, 

pi"’*’'') (150) 

The behavior of F*”''''', p*™''’'*, and their sum (F*'’'^’* + F*“'''‘*)as 
a function of internuclear distance is illustrated in Fig, 8, This shows 
that while and each approach some non-zero limit as the 

atoms separate, their sum, correctly goes to zero. The 

behavior of the components of the relaxation energies, Eqs. (144a) -(144c), 
as a function of internuclear distance is shown in Fig. 9. Qualitatively 
these correlated relaxation curves resemble those in Fig. 4 except that 
it is apparent that they fall off to zero much more quickly than their 
first-order counterparts in Fig. 4. We note that again the kinetic energy 
and the electron-nuclear attraction energies yield the largest contributions 
to the total relaxation energy. The linear relaxation in the electron- 
electron repulsion energy is shown in Fig. 9 to remain comparatively slowly 
varying with internuclear distance, the same as was found earlier for its 
first-order counterpart in Fig. 4. This indicates that this energy con- 
tribution has little or no role in determining the equilibrium bond length 
in this case. The non-linear effect, 6 is found to be negligible in 
comparison and is therefore omitted from Fig. 9. 


¥ s { 
^ ^ I 


«> 

0.05 

s. 

o 

^ 0.00 


-0.05 


- 0.10 


Inttrnucleor Distance (bohr) 

FIG. 8. The term causing the incorrect long-range behavior in the ASED 
energy through first-order, its count4rpart in the 

t(corr)^ thoi r sum as a function of 


correlation energy, F 


1 + 

internuclear distance for the Hg molecule in the X Ig state, 






Int«rnuclear Distanct ( bohr ) 


FIG, 9, The relaxation In the various components of the correlated elec 
tronic energy as a function of internuclear distance, for the H 
molecule in the state. 
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Figure 10 shows the effects on the ASED potential energy curve of 
the inclusion of correlation corrections. The two reference potential 
curves in this figure are the previously obtained (see Fig. 7) partial 
first-order ASED energy 

‘4sed’ - + <’51) 

and the full Cl potential energy curve, latter re- 
presents the best possible result attainable given the present particular 
choice of basis set. The first-order ASED energy with corrected 

by is shown in Fig. 10 to be a considerable improvement 

of the potential curve at small internuclear distances up to about 
the equilibrium distance. At larger internuclear distances, however, 
the fall-off to zero in this curve is much too slow again, and indeed 
has a peculiar shape. In contrast, the partial correlated ASED energy 
(the analog to Eq. (151)), 

I / \ nr 

(152) 


ApCI c( 0 '*‘l) _ c(corr) ^ c 4 . 

^^ASED ^ ^ ■ ae ^ 


exhibits the correct fall-off to zero at larger internuclear distances, 
but fails to attain good agreement with the full Cl curve at smaller 
internuclear distances. The two correlation effects at play here are the 
effect at smaller internuclear distances of adding and the effect 

at larger internuclear distances of the improved fall-off rate of the 
correlated relaxation energies (see Fig. 9). When both these factors 
are combined one obtains AE^^^p. This is shown in Fig. 10 to be in 
excellent agreement with the full Cl potential energy curve.^^^^ The 
small discrepancies that remain are due to the neglect of effects on 


Energy ( 10”' hartree) 


^ M 


-ASED 


.-(0 + 1) ^ _(Corr.) 

■leased + ^ 


,Ai-ci c(o-M) p(Corr.) 
'AEaseD-^ -F 


.pCO+l) 

^^ASED 


Internuclear Distance (bohr) 


FIG. 10. Comparison of ASED potential energy curves resulting from various 

levels of inclusion of electron correlation effects, and the full 

f 50 1 1 + 

Cl potential curve' S for the molecule in the X 2g state. 





54 


the matrix elements that are non-linear in the relv^iAt ‘on of the one- 
electron density matrix. The energy associated with the non-linear 
effects is 




(154) 


In addition to the contribution through Eq. (153), the non-linear re- 
laxation effects modify the solutions to the Cl secular equation, 

Eqs. (147) and (148), and through this they also affect the linear 
relaxation terms and In the present case, however, we find 

these non-linear relaxation contributions to be negligible. 

G. The Heg Molecule 

The LiH and H 2 molecules are examples, respectively, of the most 
favorable cases and the most complicated cases for the application of the 
ASED energy decomposition. There exist also cases of intermediate com- 
plication in which the ASED first-order energy does approach the separated 
atoms limit correctly, but for which the exchange energy appears in the 


least tractable form because the molecular orbitals are symmetry functions. 
An example illustrating this intermediate case is found in the Heg molecule 
in the state. The ASED zero-order wavefunction is 


= A^lOg la^agap' 


(155) 
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where the ASED occupied molecular orbitals approach the limits 


’'s) 

■ '*b) 

as the He atoms separate. The two helium atom orbitals are labeled 
Is^ and Isg respectively. The exchange energy in the He^ molecule case 
is 


E*'* = - K, , - 2K, , - K, , 

xcA 


(157) 


The linear relaxation of the exchange energy in this form is 


411 f (^l„o"K,„j)6<')(r.r')dc (158) 

r'=r 

rv 

where the limit orbitals, Eq. (156), are denoted by 10g° and lcr^^°. 

The exchange operators appearing in Eq. (158) are yet fully molecular- 
type operators. The practical advantages that derive from being able 
to express the relaxation energy as the expectation value of atomic- 
like operators are non-existent in the present form. However, one can 
reexpress the ASED energy in terms of an equivalent set of orbitals. A 
convenient equivalent set of orbitals is defined as 


IOa = (l^g + 

(^'"g ■ ''^u) 


(159b) 
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These orbitals have the property that, as the atoms separate, they 
approach the atomic orbitals themselves, 


'“A '®A 

(160a) 

loij— •U b 

(160b) 

The exchange energy in terms of these orbitals is 


4ch " ■ ‘^loA,la^ " 

and the linear relaxation in the exchange energy is 

(161) 


(162) 


We note that, by the choice of this equivalent set of orbitals in terms 


of which to construct the exchange energy, one can once again express 
the ’linear relaxation in the exchange energy as the expectation value 
of atomic-like exchange operators. The relaxation in the one-electron 
density matrix in this case is 


where 



and 

p'“'’(x-!:') = K’=a-'=a)*205b-'^b) 


(163c) 
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The total linear relaxation energy is again the expectation value of 
a Fock operator 

r'«r 

rsi# 

where, as a result of the transforniation to the equivalent orbitals, 
the Fock operator contains superpositions of atomic effective potentials, 


A. 

F 


- 1/2 + 





(165) 


D. General properties 

We conclude with some observations about the general properties 
of the ASED energy decomposition. The calculation of the ASED molecular 
orbitals, Eqs. (48)-(58), appears to require only the evaluation of one- 
electron type of integrals, while avoiding the calculation of the 
numerous two-electron multi-center type of integrals. The calculation 
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of the ASED molecular energy through first-order, Eq, (91), also appears 
to offer similar practical advantages over the calculation by the tra- 
ditional approaches of the molecular energy. The E^ component of the 
ASED energy Is a pair-wise additive term. Each pair contribution requires 
at most the calculation of two-center two-electron Integrals. The tech- 
niques dealing with this particular type of calculation have been ex- 
tensively developed and refined.^^^""^®^ Moreover, these two-center 
energies in E^ can be calculated prior to their use in a polyatomic 
molecule calculation. In many cases, the same practical advantages exist 

in the calculation of by partitioning it into E°^^ and the relaxation 
n 1 2 M 1 

terms and 6 This advantage may only exist, however, after 

reeXpressing E^^J^ in an equivalent set of orbitals. The linear relaxation 
energies, Eqs. (92) and (93), are just expectation values of Fock-operators 
containing superpositions of short-ranged atomic-like potentials. This 
lends itself to the use of matrix-representations of these atomic-like 
potentials in an atomic basis set, in a manner analogous to the procedures 
discussed in solving for the molecular orbitals. The practical problem 
of evaluating the linear relaxation energy thus reduces to the calculation 
of the projection of the relaxation in the one-electron density-matrices 
onto the various atomic orbital spaces. Provided the non-linear relaxation 
effects remain minor, the calculation of the molecular energy via the 
ASED approach has the prospect of avoiding the most time-consuming 
practical aspects of the traditional approaches. Complications, however, 
remain yet to be fully resolved in those cases when a limited Cl is re- 
quired because the ASED energy, through first-order, does not approach 
the separated-atoms limit correctly. On the other hand, it is apparent 
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that each of the Cl hamiltonlan matrix elements may themselves be decom- 
posed as we have done for the zero and first order ASED energies. The 
latter are in fact just a prototype of a diagonal hamiltonlan matrix 
element. Such a decomposition of the Cl matrix-elements, aided possibly 
by first formulating these elements in terms of equivalent orbitals, 
has prospects of resolving the problem with the correct long-range behavior 
of the ASED energy. Moreover, such an ASED decomposition of the Cl 
hamiltonian matrix has some interesting simplifying aspects in its own 
right for use in Cl calculations. The tailoring of the molecular orbitals 
to lead to a convenient decomposition of the total molecular energy is 
well known from its use in the PCIiO method. The use of the 
properties of the solutions to the atomic and diatomic subproblems for 
simplifying the calculation of the polyatomic molecular energy as in the 
present ASED method is also commonly used in the well known approach of 
atoms-and-diatomics-in-molecules.^^^*'^^^ Finally, we note that the ASED 
energy decomposition may be applied without further change to just the 
valence electrons by using, for example, the method of effective core 
potentials^^^^ to remove the core electrons from explicit consideration. 
Relativistic effects on the valence electrons may also be incorporated 
into the ASED energy decomposition without further change via the effective 
core potentials. 
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V. CONCLUSION 

We have presented a decomposition of the molecular energy. This 

particular decomposition is motivated by Anderson's ASED model of chemical 

bonding. The energy through first-order is found to adopt the 

form of the classical electrostatic energy, of a zero-order exchange energy, 

and of various relaxation energies. The latter give the effect on the 

energy of the relaxation in the one-electron density-matrix relative to 

the superposition of atomic one-electron density-matrices. We have illus- 

trated three different possible cases using the LiH X 4> , Hp X . and 
1 + a 

Heg ^ molecules as examples. The ASED approach presented here has 

the advantage of dividing the energy into physically transparent group of 

terms. Moreover, the form of the terms in the ASED energy lends itself to 

the fullest utilization of properties of the atomic solution in order to 

simplify the calculation of the molecular energy. The result is that at 

most two-center two-electron type of integrals are needed for a calculation 

if non-linear relaxation effects remain of minor importance. Thus* the 

ASED approach presented here has prospects of reducing even by an order of 

magnitude the effort required in the calculation of the molecular energy 

f 1 2 1 

by the conventional approaches.' * ' 
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FIGURE CAPTIONS 

FIG. 1. The electrostatic energy Ep, the exchange energy g, and the 
total energy change AE as a function of internuclear distance 
for the H 2 ^ molecule in the X^Eg state. 

FIG. 2. Comparison of various ASED potential energy curves for the Hg 

2 4 * 

molecule in the X Eg state with the minimum basis set variational 
energy^^^^ and the exact potential energy curve. 

FIG. 3. The zero-order electronic energy aE^^, the classical electrostatic 
energy E^, and the net ASED zero-order potential energy, 
as a function of internuclear distance for the Hg molecule in tne 
X^Sg state. The exact potential energy curve^^^^ is also shown 
for purposes of comparison. 

FIG. 4, The relaxation in the various components of the electronic energy 
as a function of internuclear distance, for the Hg molecule in the 
X^Eg state. 

FIG. 5. Decomposition of the partial ASED energy through first-order, 

AE^ggp^ - into the classical electrostatic energy, E^p , 


( 0 + 1 ) 


c>£ 


and the relaxation energy, 6E^ S for the Hg molecule in the 
X^'g"*^ state. 

FIG. 6. Decomposition of the partial ASED energy through first-order, 
AE^ggg^ - into an effective electrostatic energy, 

^ relaxation in the kinetic energy, 

for the Hg molecule in the X^2g^ state. 
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FIGURE CAPTI ONS 
(Conti nuecf) 


FIG. 7. Comparison of the partial ASED energy through first-order, 

AE^gtp^ - with the Hartree-Fock^^^^ and the exact^^^^ 

potential energy curves, for the H 2 molecule in the X Sg state. 
FIG. 8, The term causing the incorrect long-range behavior in the ASED 
energy through first-order, its count4rpart in the 

correlation energy, and their sum as a function of 

internuclear distance for the Hg molecule in the X Eg state. 

FIG. 9. The relaxation in the various components of the correlated elec- 
tronic energy as a function of internuclear distance, for the Hg 
1 + 

molecule in the X 2g state. 

FIG. 10. Comparison of ASED potential energy curves resulting from various 

levels of inclusion of electron cor -elation effects, and the full 

Cl potential curve^^^^, for the molecule in the X^ 2 ^ state. 
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